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ABSTRACT 

This paper presents results derived from analysis of new spectroscopic and 
photometric observations of the chromospherically active binary system CF 
Tuc. New high- resolution spectra, taken at the Mt. John University Ob- 
servatory in 2007, were analyzed using two methods: cross-correlation and 
Fourier-based disentangling. As a result, new radial velocity curves of both 
components were obtained. The resulting orbital elements of CF Tuc are: 
aisim=0.0254±0.0001 AU, a2sim=0.0228± 0.0001 AU, A/isim=0.902± 0.005 
Mq, and M2sim=1.008 ± 0.006 Mq. The cooler component of the system 
shows Ha and Call H & K emissions. Using simultaneous spectroscopic and 
photometric observations, an anti-correlation between the Ha emission and 
the BV light curve maculation effects was found. This behaviour indicates 
a close spatial association between photospheric and chromospheric active 
regions. Our spectroscopic data and recent BV light curves were solved si- 
multaneously using the Wilson-Dcvinncy code. A dark spot on the surface of 
the cooler component was assumed to explain large asymmetries observed in 
the light curves. The following absolute parameters of the components were 
determined: A/i=l.ll ± 0.01 Mq, M2=1.23 ± 0.01 Mq, i?i = 1.63 ± 0.02 Rq, 
i?2=3.60±0.02 Rq, Li=3.32±0.51 and 2.2=3.91 ±0.84 Lq. The primary 
component has an age of about 5 Gyr and is approaching its Main Sequence 
terminal age. The distance to CF Tuc was calculated to be 89 ± 6 pc from 
the dynamic parallax, neglecting interstellar absorption, in agreement with 
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the Hipparcos value. The orbital period of the system was studied using the 
O — C analysis. The O ~ C diagram could be interpreted in terms of either 
two abrupt changes or a quasi-sinusoidal form superimposed on a downward 
parabola. These variations are discussed by reference to the combined effect of 
mass transfer and mass loss, the Applcgate mechanism and also a light-time 
effect due to the existence of a third body in the system. 

Key words: binaries: eclipsing - stars: fundamental parameters - technique: 
spectroscopy - technique: photometry - stars: individual (CF Tuc) 



1 INTRODUCTION 



Recently, the study of chromo spherically active binaries (hereafter CAB) ha s become more 
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(119851 ) and later series of papers). CAB stars are usually either detached 
or semi-detached systems and their components include late spectral types of F, G or K, 
with luminosity classes V or IV. They show large asymmetries in their light curves, that 
are usually explained in terms of cool/dark spot(s) models. The chromospheric activity is 
associated with Call H and K or/and Ha emission lines, ultraviolet excess, soft X-rays 
and radio emission. There is a wealth of associated phenomena - including starspots, plages, 
flares, non-radiatively heated outer atmospheres, activity cycles, and deceleration of rotation 
rates. 

The study of such active binary stars provides understanding of the origin, evolution, 
and effects of magnetic fields in cool stars. In this context, we have chosen CF Tuc as the 
target of the present work. A review of this binary is presented below, new spectroscopic 
observations and their reductions are described in Section 2. The procedure used to obtain 
radial velocities and the orbital solution are outlined in Section 3, while rotation velocities 
of the components are examined in Section 4. Magnetic activity indicators. Ha and Call H 
& K emission lines, are examined in Section 5, while simultaneous solutions of the BV light 
and radial velocity curves are given in Section 6. The orbital period analysis is reviewed in 
Section 6, and a discussion of the new results is given in the last section. 

CF Tuc (HD 5303 = HIP 4157, V = 7.60 mag) is a well known southern hemisphere 
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CAB star: CF Tuc 
CAB. In the third edition of the CAB catalogue, recently published by 



EkeretaL 



(120081 ). 

the following physical properties are given: (i) It is a RS C Vn-type SB2 system with the 
spectral types G0V/IV+K4IV/V (jCutispoto and Letd (119971 )). (ii) Photometric light curves 



have remarkable asymmetries, especially visible well out of ec 



asymmetries are variable with time (e.g. iBudding and ZeilikI (119951 ): 



i pses. The magnitudes of the se 
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(iii) Radial velocities of the components were obtained, and then spectroscopic orbital ele 
ments and absolute parameters of the components were de rived: Mi =1-06 M^^ , M-^=1.21 M(^,_ 
R^=l m Rq, /22=3.32 Rq, Ti=5980 K and T2=4210 K Jcolher et all (ll98lUcoates et al 



( 1983 )) . (iv) The orbital p e riod is 2.797672 d a ys, however i t is ch anging (IThompson et al 
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i s clos e to 70°, therefore eclipses are partial (e.g. iBudding and ZeilikI (119951 ) 



(l2007al)). (v) Orb i tal inclination 



Anders et al. 



(119991 )). (vi) The projected rotational velocities of t he component s are a round 25 - 30 km/s 



for the primary, 52 - 75 km/s for the secondary (e. 
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Coates et a 



Randich et al. 
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(I2OO0I) ). (vii) The spectra 



I993I)) . (viii) The system 



Hearnshaw and Oliver! (119771 )). The Ha line ap- 



1983)). 



(ix) The system is also 



Pravdo et al 



(119961 )). Some large 



(119961 )). Other references can be found 



a known radio and X-ray sour ce (IBudding et al 
flare events were reported (e.g. iKiirster and Schi 
in the CAB catalogue. 

Particular aims of the present study are: (i) to investigate the orbital period changes 
of CF Tuc with updated data, and (ii) to derive more accurate absolute parameters of the 
system using our new, high resolution spectroscopic data. 



2 SPECTROSCOPIC OBSERVATIONS 

The spectroscopic observations of CF Tuc were made at Mt John University Observatory 
(hereafter MJUO) in New Zealand in the summer season of 2007. A Im McLennan tele- 
scope and HERCULES (High Efficiency and Resolution Canterbury University Large Echelle 
Spectrograph) spectrograph were used. The spectrograph has 100 echelle orders, which cover 
wavelength range between 380 nm and 900 nm, and is located in a vacuum tank on stable 
ground in a thermally isolated room and attached to the telescope by fibre cables. There 
are two different resolving powers: R=41000 and R= 70000. In the observations the former 
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Table 1. Journal of spectroscopic observations of CF Tuc. Signal-to-noise (S/N) ratio refers to the continuum near 5800 A. 



No 


Frame 


HJD 


S/N 


Exp. Time 






(+2450000) 




(s) 


1 


w4ooUU2z 


54d49.9SD52 


60 


1509 


2 


W4350U3U 


54350.03238 


74 


1223 


3 


W4350U53 


54350.23328 


78 


1322 


4 


w4351Ul2 


54350.89945 


85 


1144 


5 


w4o51U14 


54350.9153d 


80 


1021 


6 


w4352UlU 


54351.90655 


50 


1500 


7 


w4oo2Ul9 


54351.96099 


78 


1500 


8 


w4352U22 


r'jinr'"! noion 

54351.98189 


85 


1200 


9 


/I o r* o nnn 

w4353U02 


54352.86852 


85 


1500 


10 


w4354017 


54353.93546 


80 


1300 


11 


W4354U19 


54353.95324 


85 


1300 


12 


W4354037 


54354.10935 


95 


1300 


13 


W4354039 


54354.12710 


88 


1300 


14 


W4356003 


54355.85715 


55 


1600 


15 


W4356005 


54355.87984 


60 


1700 


16 


W4356018 


54355.99775 


70 


1500 


17 


W4356028 


54356.07125 


80 


1500 


18 


W4356030 


54356.09149 


80 


1500 


19 


W4363011 


54362.85380 


98 


1600 


20 


W4363027 


54362.93923 


80 


1900 


21 


W4363029 


54362.96309 


85 


1750 


22 


W4377014 


54376.89666 


94 


2100 


23 


W4377022 


54376.96563 


78 


2200 


24 


W4378016 


54377.94931 


82 


2100 



was used as bein g bette r adjusted to the mean seeing value [9 ~ 3.5") at MJUO given by 



Hearnshaw et al. 



mm . 

24 spectra were obtained over 10 nights between September and October 2007. The jour- 
nal of observations is shown in Table [H The exposure time was chosen between 1100 and 
1800 seconds, depending on the weather conditions. During the observations, comparison 
spectra of a thorium-argon arc lamp were taken before and after each stellar image. A set 
of white lamp images was also taken as flat field images. Two lAU radial velocity standard 
stars (HD 36079 and HD 693) were observed. The bright, non-active and slowly rotating 
standard star HD 36079 (G5II, ^^=-13. 6 km/s) was chosen for the RV measurements of 
the components of CF Tuc. Hercules Reduction Software Package (HRSP, ver. 3: Skuljan 
and Wright 2007) was used for reductions of all observations. This procedure takes into 
account other sources of apparent motion (such as the Earth's rotations and orbital revo- 
lution), wavelength calibration, removal of spurious pixels (such as those struck by cosmic 
rays), correction for the inherent pixel-to-pixel response variation (fiat-fielding), continuum 
normalization, etc., and produces a target spectrum (wavelength versus flux) as the result. 
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Table 2. Spectral orders and stellar lines used in RVs measurements of CF Tuc. 



Order No Wavelength Dominant 

Interval (A) Spectral Lines 

85 6640-6740 Sill (6660.52 A), Sill (6665.0 A) 

Fel (6677.989 A), Sill (6717.04 A) 
88 6430-6440 Fel (6430.844 A), Cal (6439.075 A), Cal (6449.81 A) 

Call (6456.87 A), Cal (6462. 57A) 
97 5820-5900 Nal D2 (5892 A), Nal Dl (5898A) 

110 5151-5188 Mgl (5174.13 A), Fel (5168.897A), Fell (5169.03 A) 

Fell (5171.595 A), Mgl (5172.6843 A), Mgl (5174.13 A) 



3 RADIAL VELOCITIES AND THE ORBITAL SOLUTION 

Measurements of radial velocities were done by two methods: The cross-correlation technique 
(CCT) and the Fourier disentangling technique (KOREL). CCT was used as the first step 
to estimate the orbital parameters, as KOREL may not produce a unique solution in some 
complex cases. In this study, the real orbit improvement comes from usage of KOREL and the 
number of spectral orders is not so relevant to accuracy. In general, if all the spectral orders 
were of the same quality, the precision of a measurement would increase with the square 
root of the number of orders. If the best spectral orders are chosen, others will produce 
systematic errors, and improvement will not be scaled as the square root of the number of 
orders. The main issue is the role of noise and better line resolution and identification in the 
selected orders. In the case of CF Tuc, since CF Tuc is rather faint for the used observational 
instruments and its components are close to GO V/IV + K4 V/IV, the observed spectra of 
the system include too many blended metallic absorption lines. It is very difficult to identify 
real/true lines and to resolve them into two components. All spectral regions were examined, 
and, as a result, four spectral orders were selected, for which the lines of both components 
could be clearly detected, and therefore suitable for disentangling. The information about 
these four spectral orders is given in Table El 



For the CCT, the FXCOR tas k in the radial velocity package of IRAF (ITonry and Davis 



(I1979I ): iPopper and Jeongl (jl994l )) was used. FXCOR calculates the velocity Doppler shift 
between two spectra (of the variable and comparison stars) by fitting the correlation with 
a user-selected function. In the present study, the Gaussian function was adopted as the 
best-fitting one. The spectra of HD 36079 were used as a template for deriving RVs of the 
components. In order to obtain orbital parameters from the radial velocity data derived from 
the CCT, the ELEMDR77 program, developed by T. Pribulla (2008, private communica- 
tion), was used. 
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6430 6440 6450 6460 6470 

Wavelenght(A) 

Figure 1. KOREL solution of cchcUc order 88 as a sample. Upper panel shows the KOREL fits while the lower panel displays 
the disentangled spectra of each component. 



For the Four ier disentangling technique, the KOREL code (developed by lHadraval (119951 ). 



Hadraval (119971 )) was used. In the first step with this procedure, KOREL requires input 



parameters within realistic bounds. The initial values of parameters were taken from the 
solution of RV curves obtained with the OCT method. Four spectral orders, given in Table 
[21 were analyzed simultaneously. After several iterations, the KOREL code gave a value 
close to for the eccentricity e within its uncertainties. For this reason, we assumed a 
circular orbit for the system. Additionally, during the fitting, the orbital period Porb of the 
system was fixed to be 2.7975004 days (see Section 5). The velocity amplitudes Ki and K2 
of the components and the conjunction time Tq were the adjusted parameters. The best 
fitting orbital elements are given in Table [3], and the best fits to the composite spectra and 
disentangled spectra are shown in Fig. [T]for the echelle order 88 as a sample. 

The KOREL code can not derive the systemic velocity of the binary star, however, KO- 
REL retains the systemic velocity in the disentangled spectra of the components. Therefore, 
the systemic velocity was adopted by averaging systemic velocities obtained from the 
CCT method and the ELEMDR77 program. The adopted systemic velocity was added to 
the RVs measured by KOREL to yield the final RVs of the components, which are given in 
Table m In this table, phase values of observed time of RVs in the second column were calcu- 
lated using the linear ephemeris given in Eq. 1. 0-C values in the fourth and sixth columns 
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Table 3. Spectroscopic orbital parameters of CF Tuc 



Parameter 


Value 


Porb (days) 


2.7975004 (fixed) 


To (HJD) 


54327.0583±0.0012 


y^(km/s) 


9.58±0.14 


q 


1.117±0.009 


Xi(km/s) 


98.92±0.24 


A'2(km/s) 


88.55±0.24 


aisinj(AU) 


0.0254±0.0001 


a2sini(AU) 


0.0228±0.0001 


Misini(M0) 


0.902±0.005 


M2sini{MQ) 


1.008±0.006 



Table 4. RV measurements, with O-C values from theoretical fit, of components of CF Tuc. 



Time 


Phase 


RVi 


O-Ci 


RV2 


O-C2 


HJD 


</- 


(km s^-"-) 


(km s^-"-) 


(km s-i) 


(km s-i) 


2454377.9493 


0.192 


-82.4 


1.0 


91.8 


-0.5 


2454349.9865 


0.196 


-84.6 


-0.4 


92.3 


-0.8 


2454350.0324 


0.212 


-87.2 


-0.2 


96.2 


0.8 


2454352.8686 


0.226 


-87.9 


0.7 


96.8 


0.0 


2454350.2333 


0.284 


-87.7 


-0.7 


95.5 


0.1 


2454355.8572 


0.295 


-85.3 


0.0 


94.4 


0.6 


2454355.8798 


0.303 


-84.4 


-0.7 


92.6 


0.1 


2454355.9978 


0.345 


-72.5 


-0.7 


87.0 


4.7 


2454356.0713 


0.371 


-62.1 


-0.6 


73.7 


0.3 


2454356.0915 


0.378 


-58.9 


-0.5 


70.8 


0.1 


2454350.8994 


0.522 


24.9 


0.5 






2454350.9154 


0.528 


26.7 


-1.3 






2454353.9355 


0.608 


71.1 


-1.7 


-45.4 


0.9 


2454353.9532 


0.614 


75.7 


0.1 


-48.2 


0.6 


2454354.1093 


0.670 


95.8 


-1.3 


-67.5 


0.6 


2454354.1271 


0.676 


96.0 


-2.8 


-71.8 


-2.2 


2454362.8538 


0.796 


104.3 


-0.1 


-75.1 


0.3 


2454376.8967 


0.815 


100.2 


0.0 


-70.4 


1.4 


2454362.9392 


0.826 


97.3 


0.0 


-68.8 


0.3 


2454362.9631 


0.835 


94.8 


0.3 


-66.6 


0.1 


2454376.9656 


0.840 


93.0 


0.2 


-65.0 


0.2 


2454351.9065 


0.882 


76.3 


0.4 


-50.0 


-0.1 


2454351.961 


0.902 


66.9 


0.9 


-41.7 


-0.5 


2454351.9819 


0.909 


63.1 


0.7 


-38.2 


-0.4 



represent the residuals between observed and theoretical RVs obtained from simultaneous 
solution of the BV light and RV curves, described in Section 5. 



4 ROTATIONAL VELOCITIES 



The program PROF (IBudding and ZeilikI (119951 )). which follows an ILOT type curve- fitting 



procedure, was used to determine rotational velocities. 



tional broadenings, as discussed by 



Budding and Zeilik 



'RO 



convolves Gaussian and rota- 



19951 ) , and computes the line profile 



as a function, basically, of the following parameters: the continuum intensity Jc, the relative 
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5889 5890 5891 5892 5893 5894 

Wavelength 



5889 5890 5891 5892 5893 




5889 5890 5891 5892 5893 

Wavelength 



Figure 2. Results of profile fitting to Na D2 5892 lines at phase 0.371 (upper panel) and phase 0.826 (lower panel). 



Table 5. Profile fitting parameters at two different orbital phases for the Na D2 5892 feature. 





phase 


0.371 


phase 


0.826 


Parameter 


Primary 


Secondary 


Primary 


Secondary 


Ic 


1.024±0.015 


0.990±0.010 


0.927±0.012 


1.031±0.008 


Id 


0.216±0.018 


0.188±0.010 


0.192±0.019 


0.175±0.008 


A„(A) 


5888.740±0.043 


5891.452±0.094 


5891.852±0.043 


5888.641±0.102 


r(A) 


0.498±0.053 


1.142±0.108 


0.528±0.055 


1.216±0.088 


s(A) 


0.469±0.107 


0.745±0.074 


0.300±0.060 


1.053±0.055 


r(km/s) 


25±3 


58±6 


27±3 


62±5 


s(km/s) 


24±6 


38±4 


15±3 


54±3 




0.01 


0.01 


0.01 


0.01 




1.034 


1.040 


1.024 


1.048 



depth Id at mean wavelength Am, the rotational broadening parameter r, Gaussian broad- 
ening paramete r s of a given line, a nd the limb darkening coefficient u. A simila r procedure 



was followed by 



OlahetaL 



fll992h . 



Olahetal 



(119981) and 



Budding et al 



(120091). 



We considered the Na D2 line profiles for CF Tuc in the Hercules spectral order 97 and 
fitted the selected line profiles at various orbital phases using PROF. Typical results of the 
profile fitting at phases of 0.371 and 0.826 are shown in Fig. [2] and given in Table [51 

According to the value of r (the rotational broadening parameter) in Table [5l the 
projected rotational velocities of the primary and secondary components are 26±3 and 
60±5 kms~^, respectively. Using absolute parameters of components from Table [H and 
frotsin2=27ri?sim/Prot (assuming synchronous rotation, Prot=Porb, and irot=iorb)^ theoretical 
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rotational velocities were be found as 28±1 and 62±2 kms~^ for the primary and secondary 
components, respectively. Therefore, we found that within the error limits, synchronous 
rotation for both components can be reliably accepted. 

The Gaussian broadening parameter s (presented in Table [5]) ranges from 24 to 15 km/s 
for the primary component and from 38 to 54 km/ s for the secondary component at phases of 
0.371 and 0.826, respectively. This parameter generally relates thermal broadening and other 
broadening factors (micro and/or macro turbulence, etc.) rather then rotational broadening. 
If the temperature of 4300 K is taken, where the Na D2 lines are formed in a subgiant 
atmosphere, thermal velocities are found to be ~2 km/s. Therefore, for this value, the 
thermal broadening is not significant, and the broadening parameter must include some other 
motions, possibly related to turbulence effects or magnetic activity. Since our spectroscopic 
observations and Innis's CCD BV observations were made almost simultaneously, we can 
compare these spectroscopic results with the spot model derived from the BV light curve 
analyses. For instance, when the maculation is close at phase 0.821, the s parameter is 
greatest and falls to about 38 km/s half a cycle later (at phase 0.371). Therefore, we can 
deduce that the larger value of s of the secondary component could be associated with the 
relatively stronger effects of surface activity of this star. 



5 MAGNETIC ACTIVITY INDICATORS 

The Ha and Call H & K emission lines are very important indicators of magnetic activity - 
in other words, chromospheric activity. Generally, the more active stars show these emission 
lines always above the continuum (e.g. UX Ari, II Peg, AR Psc, V711 Tau and XX Tri). 
Except for being sensitive to the chromospheric activity, these emission lines are also a 
good diagnostic of inter-components matter in the form of gas streams, transient or classical 
accretion disks and rings in mass-transferring binaries i.e. the Algol t ype, in which the cooler 



star fi lls its Roche lobe and transfers mass to the hot companion (e.g. iRichards and Albright 



mm ). 

The Ha and Call H & K observed line profiles of CF Tuc are shown in Figs. [3], [Hand [51 In 
these spectra, the Ha profiles show absorption and emission features, while the Call H & K 
are in emission at all orbital phases. All absorption and emission features were red and blue- 
shifted depending on the orbital phases. Since the cooler component of the system is an active 
star, it is assumed that these emission features are related to this component. To confirm 
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Figure 3. Observed Ha line profiles for CF Tuc in 2007 at first half (left panel) and second half (right panel) of the orbital 
period. The vertical lines represent the laboratory wravelcngth (A()=6562.82 A) of the Ha line. 



this, we calculated the RV values of the Ha and Call K emission features and plotted them 
in Fig. [6] with RVs of both components of the binary system. We took into account only the 
K line (3933.66 A) here, because of likely contamination from He centered only 1.5 A away 
from the H line (3968.47 A). As it can be seen in Fig. [6], although not strictly sinusoidal, the 
RVs of the Ha emission feature follows the orbital motion of the cooler component but with 
a larger amplitude of about 200 kms"^. This rather larger amplitude indicates that the Ha 
emission feature could originate in a gas cloud in the form of a chromospheric prominence 
from the cooler component. The RVs of the Call K emission line closely follows that of the 
cooler component (see Fig. [6]), indicating its origin to be in the chromospheric layers of that 
star. 

If we assume that the Ha emitting region is at rest in the rotating reference frame of 
the system and that it lies between th e barycentre of the s ystem and the cooler component. 



we can estimate its location following 



Marino et al. 



200ll ). who modelled inter-components 
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Figure 4. Sample of observed Call H & K spectra for CF Tuc at different orbital pliases. 



matter using Ha lines for HR 7428. The HR 7428 star is an RS CVn-type binary, composed 
of a bright K giant and an A- type Main Sequence dwarf. It is a well detached system like 
our target star, CF Tuc. In this model, the projected distance of the emitting region can be 
found from aer=cic{Ker/ Kc), where Cc is the semi-major axis of the cooler star orbit, and 
Kc and K^^r are the semi-amplitudes of the RV variations for the cooler star and emitting 
region, respectively. Since our orbital solution gives Kc=89 kms~^ and ac=3.6xl0® km, we 
find ae,.=8.2xl0^ km - in other words, the Ha emitting region should be located at a distance 
of about 2 Rc from the surface of the cooler component. We estimate almost the same value 
for the location of the emitting region using Ver,rotsini=2naer sini/Prot (with the assumption 
of synchronous rotation). 

The simultaneous photometric and spectroscopic observations of CF Tuc offer the possi- 
bility of studying the photospheric (spots) and chromospheric (plages and/or prominences) 
active regions of the cooler component. With this aim, we measured the equivalent widths 
[EWs) of emission features observed in the Ha and Call K spectra by means of multiple 
Gaussian fits, and plotted them versus orbital phase in Fig. [71 In Section 6, we present a 
large cool photospheric spot on the cooler component to explain observed BV light curve 
asymmetries. To concentrate on effects due to the spot only, we subtracted the eclipse and 
proximity effects from the observed data using the unspotted (immaculate) light curve pa- 
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3925 



3930 



3935 



3940 



Figure 5. Observed Call K line profiles for CF Tuc in 2007 at first half (left panel) and second half (right panel) of the orbital 
period. The vertical lines represent the laboratory wavelength (Ao=3933.66 A) of the Call K line. 




Figure 6. Variations of RVs of emission features observed in Ha and Call K. The continuous lines are the RV solutions for 
the primary (hotter) and secondary (cooler) components, according to parameters given in Table 6. 
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Figure 7. Maculation effects in BV magnitudes (upper two panels), and observed EWs of H« and Call K emission features 
(lower two panels) as a function of tlic orbital phase for CF Tuc. 



rameters given in Tableland plotted these points (which show the maculation/distortion 
wave) versus the orbital phase in the upper panels of Fig. [71 

The maximum EWna value of about 1 A is reached at phases 0.6 - 0.7, where the spot 
effect is dominant {BV distortion of lights minima); the minimum value of EWua —0.2 
A is observed between phases 0.0 and 0.5, where the spot could not be seen in the light 
curves. The Call K equivalent width does not show any clear orbital modulation (rotational 
modulation under the synchronous rotation). However, there is a possible ant i- correlation 
between BV light curve asymmetries and Ha emission, which is apparent with an almost 
similar shape of the curves. This behaviour denotes a close spatial association of photo- 
spheric and chromospheric active regions. Such anti-correlations between photospheric and 



chromospher ic diagnostics are fo u nd in some a c tive s t ars and have 



authors (e.g. 



Frasca et al. 



fl2000h . 



Frasca et al. 



fl2005h : 



Biazzo et al. 



j een ex amined by many 
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6 SIMULTANEOUS SOLUTION OF THE BV LIGHT AND RADIAL 
VELOCITY CURVES 

We analyzed the BV light curves from Innis, J. L. (2008, private communication), Hipparcos 
light curve (ESA 1997) and radial velocity curves f rom th is study using the Wilson-Devinney 
code (WD), version 1996 (jWilson and Devinneyl (jl97ll )). The BV light curves from Innis 
and our new radial velocity curves were solved simultaneously. Innis et al. observed CF 
Tuc in BV filters at Brightwater Observatory in the summer season of 2007. They used 
a short-focus, 70mm telescope and a cooled SBIG ST7E CCD camera, which gives a field 
of view of near 0.8 arcdegreexO.55 arc degree. A detailed description of the observatory 
and techniques is given in the paper by Innis et al] (2007b). They observed HD 5210 and 
HD 4644 as comparison and check stars, respectively. We calculated external uncertainties 
for all comparison minus check magnitudes and found them to be 31 and 13 mmag in B 
and V filters, respectively. For this, we used the standard deviation of the differential light 
variations of the comparison relat ive to the ch e ck sta r collecte d during the same night . A 



similar procedure was followed by 



Erdem et al. 



fl2009h and also 



Strassmeier et al. 



fll999h to 



examine the quality of long-term multicolour photometric data of several active stars. The 
observational data were not transformed into the standard BV system. It is worth noting 
that our spectroscopic observations and their photometric observations were made almost 
simultaneously. In order to calculate the phases of the CCD BV light observations of CF 
Tuc, the light elements of the system were derived by using the photoelectric primary minima 
times with E >2000 cycles (see Fig. [Qb) as; 



HJD {Mm I) = 2452452.7326(68) + 2*^.7975004(9) x E. 



(1) 



with the weighted least squares method. In the light curves which were formed using these 
phases, the primary minimum coincides with the phase 0.0 (see upper panel of Fig. [8]). 

In the WD method, some parameters could be fixed according to theoretical models. In 
the light curve modelling, the temperature of the primary component was fixed at 6100 K, fol- 
lowing Anders et al.(1999) and Budding and McLaughlin (1987). The r oot square limb dark- 



ening law was adopted, a nd th e darkening coefficients were taken from 



Diaz-Cordoves et al 



(119951) and 



Claret et al 



(119951 ) The bolometric gravity-d arkening coef ficients of the compo- 



nents were set to 0.32 for convective envelopes, following 



Lucy 



albedos were fixed to 0.5 for convective envelopes, following 



1967) 



also, the bolometric 



Rucinskil (119691 ). According to 



analysis of the rotational velocities (see Section 4), the components rotate synchronously. 
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Table 6. Final solutions of light and radial velocity curves of CAB star CF Tuc 



Parameter 


BV + 2B.V 


Hp 


a{RQ) 


11.08±0.02 




Phase shift 


-0.0017±0.0002 


-0.0003±0.0006 


(km/s) 


9.6±0.4 




i (deg) 


69.91±0.09 


69.91 


Ti (K) 


6100 


6100 


T2 (K) 


4286±19 


4286 


Qi 


7.907±0.086 


7.907 


Q2 


4.452±0.011 


4.452 


Qcorr = m2 / mi 


1.115±0.003 


1.115 


I1/I12 (B) 


0.616±0.008 




h/h2 (V) 


0.557±0.008 


0.567±0.001 


n {mean) 


0.148±0.001 


0.148 


r2{mean) 


0.325±0.001 


0.325 


Spot parameters 


Spotl co-latitude (deg) 


155±5 


121±5 


Spotl longitude (deg) 


303±4 


329±3 


Spotl radius (deg) 


47±2 


31±4 


Spotl Tgpot/Tstar 


0.746±0.038 


0.731±0.026 


Spot2 co-latitude (deg) 




24±3 


Spot2 longitude (deg) 




206±3 


Spot2 radius (deg) 




38±2 


Spot2 Tspot/Tatar 




0.785±0.028 


j:w{o - C)2 


0.03177 


0.01192 



Therefore, the rotation parameters were assumed as F/i=Fc=l. From the spectroscopic or- 
bital solution described in Section 3, the circular orbit (e=0) was adopted. 

The adjusted parameters in our computations are: the semi-major axis of orbit (a), phase 
shift, systemic velocity of the binary (V^), orbital inclination (z), surface temperature of the 
secondary component (T2), non-dimensional surface potentials of both components {Qi and 
^2), and the fractional monochromatic luminosity of the primary component (li/(li+l2)). 
The initial values of q, A and were taken from the radial velocity solution (see Section 3). 
Due to the probability of the existence of a third body, resulting from the orbital period 
analysis of the system (see Section 7), a third light contribution (l^) was also considered as 
a free parameter. However, we soon found its contribution to be negligible. The binary CF 
Tuc, as mentioned above, is a typical RS CVn type eclipsing binary and its light curves show 
distortion wave like other systems in this group. In the BV light curves, the distortion wave 
appears to be an asymmetry between the light levels of the maxima (see Fig. [HK)- Therefore, 
we had to consider a cool spot on the secondary and allowed the spot parameters to be 
adjusted. 

In order to get good starting parameters for the WD code, mean points were calculated 
from all individual observations. This provided standard deviations for B and V filters; the 
highest values were 15 mmag (the clump near phase 0.9) in B, and 11 mmag (near phase 
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0.1) in V . These mean points (69 in B, 64 in V) were used with a Monte Carlo search to find 
a fit to the BV hght curves with a mass ratio fixed at the value derived spectroscopically 
(see Table [3]). The resulting parameters were next set as starting values for simultaneous 
solution of the BV light and RV curves. Simultaneous BV light and RV curves solution 
was done for all individual points, assigning the same weight and different errors for the B 
and V filters as resulted from the calculation of mean points. Errors for the RV curves were 
assigned the same as their quality was comparable. The NOISE control parameter was set 
to 1. 
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Simultaneous convergent solutions of the BV light and RV curves were obtained by it- 
erations, until the corrections of the parameters became smaller than their corresponding 
errors. The results of the final solution are given in Table El The comparison between ob- 
served and computed light curves is shown in Fig. [8^, while that of RV curves is presented 
in Fig. 6. The three-dimensional model demonstrating the presence of a large dark spot on 
the surface of the cooler component a nd the Roche geometry o f the system (making use 



Bradstreet and Steelmaru (120021 )) are also shown in 



of the Binary Maker program, ver.3.0. 
Fig. Ep,c. 

We made additional attempts to check the stability of our solution. This was done by 
assigning larger errors for the photometric light curves (especially for the B light curve). 
Additionally, only RV curves (adjusting only the parameters relevant to the orbit) were 
solved using the WD code. It was found that the simultaneous solution of BV and RV curves 
resulted in the mass ratio q=1.115±0.003, while only the RV solution gave q=1.113±0.004. 
Also, the solutions with assigned larger errors for B data gave (within uncertainties) a similar 
value of the mass ratio. 

The Hipparcos light curve of CF Tuc is available from the Hipparcos web page and 
contains 121 points with an average observational error of 11 mmag. Before starting analysis, 
Hipparcos observations of the system were t ransformed to Johnson V ma gnitudes using 



Hp — V = 0.22 x(i? — V) calibration given by lRucinski and Duerbeckl (119971 ) 



The following ephemeris was used to phase the Hipparcos photometric data: 
HJD {Min I) = 2448502.560 + 2^^.79765 x E. (2) 

The data was weighted by using the equation Wi = l/crf, where cTj is the individual 
standard error of the data given in the Hipparcos catalogue. About ten photometric points 
were discarded due to their relatively large errors. We used the mean maxima levels at phase 
0.25 for the flux normalization of both the Hipparcos and 2007 BV filters data. 

During the iterations, only spot parameters, phase shift, and luminosity of the primary 
component were treated as free parameters; others were adopted from the simultaneous 
solution of BV and radial velocity curves. As can be seen in Fig. [8^, the Hipparcos light 
curve shows two large asymmetries, one at about phase 0.70 and another at the primary 
minimum. Therefore, the possibility of two dark spots on the secondary component was 
considered. The final results are given in Table [6] and displayed in Fig. [8^,b. 
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7 ORBITAL PERIOD ANALYSIS 

In order to investigate the orbital period variation of the system, we gathered 33 minima 
times availab l e frorn the lists compiled by Kreiner, J. M. (2008, private communication) and 



Anders et al 



(119991 ) and we added one minimum time, which was calculated from the 2007 



BV light curves observed by Innis et al. (2008, private communicati on). As a first s tep, — C 



Anders et al. 



fll999h 



values were calculated using the following light elements, given by 

HJD {Min I) = 2444219.270 + 2^^.797715 x E. (3) 

The — C values versus E values (and years) we re found and the results are shown in 
the upper panel of Fig. [9K,b. [Thompson et al.l (jl99ll ) first noted that the orbital period of 



the system changes in the form of an upward parabola and t ried to explain this variation 



in terms of a mass transfer or a mass loss from the system. 



Anders et al. 



fll999h showed 



that the orbital period change has a cyclic character and discussed the — C diagram, the 



spot-wave a mplitude arid the m ean 



mechanism. 



Innis et al. 



toom and 



ight change of the s ystem as being due to the Applegate 



Innis et al. 



(j2007al ) reported that the orbital period of 



the system did not show any change between 1995 and 2006. Therefore, we could say that 
the real nature of the period variation shows up as the data increases by time. 

In the p resent — C analysis, due to a large scatter, one spectroscopic time of minimum 
obtained by lHearnshaw and Ohverl (119771 ) was discarded and altogether 33 photometric data 
were used. The standard errors of observed minima times, as given by authors, are shown as 
error bars in Fig. [9^,b. The weights were assigned according to these errors. As the standard 
errors are given in 4, 3 and 2 decimal places, we used 10, 7 and 5 for weights, respectively. 
The observed long-term period decrease of CF Tuc from these — C diagrams could thus 
be explained as follows: 

(i) Abrupt period changes: The — C diagram in Fig. [9K was considered in terms of 
abrupt period changes. Period jumps might have occurred two times within an interval of 
about 30 years. Considering that the period has remained constant between these two jumps, 
we calculated following three linear ephemerides. The first ephemeris valid for E ^ 776 is: 

HJD {Min I) = 2444219.2615(52) + 2*^.797692(12) x E, 

for 885 ^ E ^ 1835 we have: 

HJD {Min I) = 2444219.2126(51) + 2*^.797749(4) x E, 
and finally for E ^ 2192 we have: 



(4) 



(5) 
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HJD {Mm I) = 2444219.7446(153) + 2^^.797482(5) x E. (6) 

The first abrupt period change AP/P=(-2.04±0.31)xlO"^ occurred at HJD 2446527±70; 
while the second one AP/P=(+9.51 ± 0.06) x 10"^ occurred at HJD 2449745±100. The 
time interval between these two possible abrupt period changes was taken as 3218±105 
days (or approximately 9 years). Such sudden period jumps could be caused by anisotropic 



mass ejections from one (or both) component (s) (e.g. iHuang (119631 )). Indeed, the RVs of Ha 
emission support such a high velocity ejecting gas from the active component (see Section 5), 
like a prominence in the solar atmosphere. However, as seen in Fig. [9^, the sudden period 
changes seem to have occurred in the pattern of one period increase and a subsequent 
decrease. This might be an indication of sinusoidal variations rather than abrupt period 
changes due to sudden mass ejections. 

(a) Continuous period change and the light-time travel effect: A reasonable fit to the 
O — C data is obtained by using a sinusoidal ephemeris with a quadratic term, as 

C = To + PxE + QxE^ + A,sin[^{E - T,)], (7) 

where Ag is the semi-amplitude, P, the period and T, the time of minimum. However, to 
show the parabolic — C change clearly, we recalculated — C values by using the following 
light elements; 

HJD {Mm I) = 2448922.2310 + 2"^. 797630 x E. (8) 

These new O — C2 values were then plotted against the epoch number and observation years 
in the second panel of Fig. [9b. This panel shows a simple downward parabola, where the axis 
of symmetry is parallel to the O — C-axis with the vertex at epoch number E=0. It should 
be noted that the best model fitting both O — Ci and O — C2 values yields the same values 
within uncertainties for both quadratic and sinusoidal terms. Finally, parameters of the best 
theoretical curve fitting the — C data are given in Table U\ and the best theoretical fit 
with the observational data is plotted in Fig. [9)d. 

According to the quadratic term given in Table [71 the orbital period of CF Tuc is contin- 
uously decreasing at a very rapid rate of 1.11 ±0.05 seconds per year. Here, CF Tuc appears 
to have the highest rate of period decrease among the RS CVn systems. We considered 
the combined effect of the mass loss and the mass transfer to s tudy this obs e rved p eriod 



decrease of t h e syste m and used the following equation, given by lErdem et al.l (l2007a] ) and 



Erdem et al. 



fl2007bh : 
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Figure 9. (a) O — C diagram of CF Tuc constructed with Eq. (3). It shows that the period variation of the system could be 
interpreted in terms of two abrupt period changes. The residuals were calculated separately with the three linear ephemerides 
given Eqs. (4,5,6). (6) Sinusoidal representation (solid line) is superimposed on the parabolic form (dashed line) of the O — C 
variation of CF Tuc, and the residuals from the best fit curve. The O — C\ values in the top panel were calculated by using 
the linear ephemeris of Anders et al. (1999), while those of O — C2 in the second panel were computed by Eq. (8) to show a 
symmetric parabola at the vertex with E=0. These panels show the same parabola + sinusoidal variation, (c) Systemic velocity 
variations in CF Tuc. The individual points show observed systemic velocities; the continuous curve is the radial velocity curve 
corresponding to the light-time orbit in Table [7] (d) A relative astrometric orbit of CF Tuc about a barycentre with a third 
body. The points represent Hipparcos astrometric data, while the solid curve corresponds to the astrometric solution given in 
Tabled 
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Table 7. Parameters for the CF Tuc sinusoidal O-C solution + astrometry. 



J. diL (iiii^uci 




Sinusoidal O-C solution 




To (HJD) 


2448922. 2102± 0.0021 


P (d) 


2.797641± 0.0000015 


Q (d) 


-4.9 X 10-®±2.1 X 10-9 


As (d) 


0.0363± 0.0018 


Ps (yr) 


17.87± 0.57 


Ts (HJD) 


2430387± 599 


siy(o - C)2 


0.00848 


Astrometry solution 




Pl2 (yr) 


17.87 


Tia (HJD) 


2443441 


ai2 (mas) 


72.04 


ei2 





iJl2 (deg) 


90 


ii2 (deg) 


83 ± 7 


(deg) 


144 ± 19 


AacosS (mas) 


37.4±0.2 


AS (mas) 


56.4±0.2 


AfiaCosS (mas/yr) 


-2.9±0.3 


A/is (mas/yr) 


0.7±0.3 


AH (mas) 


-0.7±0.3 




1.01 



p " V o y M - M,M, <-'^ 

where the mass AM is transferred from the mass-losing component to the gainer, the 6M is 
the amount of mass lost from the system after co-rotating with the system up to the distance 
ta (i.e. Alfven radius), and AP/P is the period change. Since CF Tuc is a detached system, 
in which the primary and secondary components are filling ~ 52% and ~ 89% of their lobes 
(see Section 6), a direct mass transfer between components is not expected. However, the 
secondary component is a magnetically active and larger star, which is not far from filling its 
Roche lobe, and then there could be weak coronal flow from this component to the primary 
component through the inner Lagrangian point. Therefore, the active component might have 
a strong stellar wind, which drives the mass loss and mass transfer in the system. The RVs 
of Ha emission support such a strong stellar wind, which reaches twice the larger distance 
than the radius of the active, subgiant component (see Section 5). If we assume that the 
transferred mass due to the wind from the secondary to the primary component and the co- 
rotating distance are 10~^^ Mq/jt and 10-R2; respectively, then Eq. (9) gives the mass loss 
rate of 6M = 3.38 x 10^'' Mq/jt for the observed period change of AP/P = —4.57 x 10~^ 
yr~^. It is worth mentioning that this mass loss rate is 10 ti mes higher than the maximum 



value of the range between 10 and 10 ^ Mq/jt given by iHilditchI (120011 ) for the mass 
losses due to winds from red-giant stars. 
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There are two plausible causes of the sinusoidal O — C variation: a light-time effect due 
to a third body in the system and a period modulation due to the magnetic activity cycle 
of one of the components. We shall investigate these suggested hypotheses in turn. 

According to Table [71 CF Tuc would have a circular orbit around the center of the mass 
of a three-body system and its period would be 17.87 ± 0.57 yr. The projected distance 
of the center of mass of the eclipsing binary to that of the three-body system would be 
6.29 ± 0.31 AU. These values lead to a large mass function of /(M3) = 0.78 ± 0.07 Mq for 
the hypothetical third body. The mass of such a third body would then range from 9.63±0.29 
Mq for Zi2=30° to 2.71 ± 0.13 Mq for Zi2=90°. Here the sum of masses was taken as Mi 
-|- M2=2.34 Mq (see Section 8). If the third body were co-planar with the eclipsing pair, 
its mass and the radius of its orbit around the center of mass of the three-body system 
would be about 2.99 ± 0.13 Mq and 4.93 AU, respectively. This value of r-^ is smaller than 
the radius of the orbit of Jupiter, however, it shows that the third body would revolve far 
beyond the outer Lagrangian points of CF Tuc, and its orbit should be stable. If we consider 
the distance of CF Tuc as 89 pc (see Table 8), the minimum projected angular separation 
between the third star and the eclipsing pair could be estimated as ~ 71 mas. 

The semi-amplitude of the radial velocity of the center of mass of the eclipsing pair, rela- 
tive to that of the three-body system, is derived to be 10.5 km/s, which is a convenient value 
for modern spectroscopic observations to resolve reliably. The theoretical variation of the 
systemic velocity of CF Tuc, caused by the orbital motion around the common barycentre, 
is illustrated in Fig. [nb. There are thre e value s of the systemic velocity o bserve d at different 
epochs: 12.1±4.7 km/s (jColher et all (|l98ll)), 0.5±1.6 km/s (IBalonal ilMzl)), and 9.6±1 
km/s (present study). Except for the first data point, the observed systemic velocities follow 
the long-term variation corresponding to the light-time orbit in Table [7] and are almost the 
same (within their standard erros) as the theoretical values. The velocity measurement by 



Collier et al. 



(Il98ll ) is affected by a relatively large error, which could be caused by the 
scatter of RV data points and the two year time span of their observations. 

The astrometric method was also use d to check the t hird body hypothe sis. A similar pro- 



cedure was applied by 



Ribas et a. 



Bakis et al 



fl2005h and 



(2002 ) forRCMa, b y 

\\1 

Budding et al.l (120091 ) for U Oph and also 



Bakis et al. 



Zasche and Woli 



(2006) for XY Leo and 5 Lib, by 
(120071 ) for yW C ep, C, Phe and HT Vir. We used the Hipparcos Intermediate Astrometric 
Data (lESAl (119971 )) for this. Hipparcos observed CF Tuc between January 1990 and January 
1993. There are 76 one-dimensional astrometric measurements corresponding to 40 differ- 
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ent epochs in the Hipparcos Intermediate Astrometric Data, which were obtained by the 
two Hipparcos data reduction consortia: FAST and NDAC These data are available from 
the Hipparcos web page. In fact, the astrometric method gives support to the third body 
hypothesis on orbital period analyses of eclipsing binaries in two ways: one is to plot an 
orbit of the eclipsing binary about the barycentre of a three-body system and the other is to 



determ ine its orbital inclination (iu)- We followed the procedure applied by iBudding et al. 



(120091 ) for our target. In this procedure, an orbital model, which is derived from the orbital 
motion of the eclipsing binary around the barycentre of a three-body system, is convolved 
with the astrometric motion (parallax and proper motion). This model has 12 independent 
parameters: ai2, ei2, uJi2, ii2, Pu, ^12 (periastron passage time), (seven for the orbital 
parameters), and a, 6, fia, n (five for the astrometric components; equatorial coordinates 
+ proper motion + parallax). Since the Hipparcos astrometric data cover only 1/6 part of 
the orbital period of the hypothetical three-body system, we could take only two parameters 
from the orbital parameters, iu and Q12, together with five adjustable astrometric param- 
eters. The final results are given in Table [7] and Fig. [9ll. The orbital inclination of CF Tuc 
in the triple system is, iu, about 83° and gives the mass of a third body as about 2.74 Mq. 
Unfortunately, since the time span of the Hipparcos observations is much shorter than the 
orbital period of the three-body system and is not the periastron passage of CF Tuc in the 
three-body orbit, the data, as shown in Fig. [9ll, cover only a small part of the orbit. 

An alte r native way of explaining the 0-C behavior would be Applegate's mechanism 
(jApplegatd (119921 )). According to this, the cyclic magnetic activity could produce orbital 
period modulations, which are observed in some eclipsing binaries, especially in RS CVn 
type systems. Magnetic activity can change the quadrupole moment of a component as the 
star goes through its activity cycle. The cyclic exchange of angular momentum between the 
inner and outer parts of the star can change both the shape and radial differential rotation of 
the star. The torque required for such transfer of angular momentum could be provided by a 
subsurface magnetic field of several kC Any change in the rotational regime of a component 
of an eclipsing binary due to magnetic activity could be reflected in the orbit, as a conse- 
quence of the spin-orbit coupling. Here, we shall use Applegate's formalization to examine 
the sinusoidal part of the orbital period variation of CF Tuc and assume that the secondary 
star could be responsible for the observed orbital period modulation. The O — C diagram of 
CF Tuc, given in Fig. [9]3, shows a modulation with a semi- amplitude of 0.0363 days and a 
modulation period of 17.87 years. This gives AP/P=3.49xlO~^ The change in the orbital 
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period is AP=8.45 s. The angular momentum transfer would be AJ=6.84xlO^'^ gcm^s"^. 
If the mass of the shell is Ms=0.1M, the moment of inertia of the shell is /s=l. 03x10^^ 
gcm^, and the variable part of the differential rotation is Af2/r2=0.026. The energy bud- 
get and RMS luminosity variation are A£'=9.08x 10^^ ergs and ALj:ja,/5=13.23 Lq. This 
model gives a mean subsurface field of 12.1 kG. The RMS luminosity variation predicted 
by this model is larger than the total luminosity of the active star. Therefore, a model 
with M.s=0.1M a nd VLcir=^^ cannot explain the orbital period change observed in CF Tuc. 



Applegatd (119921 ) has also calculated a similar result for RS CVn itself. However, he could 



obtain a reasonable result for RS CVn using the following two approaches: one is that the 
active component has solid body rotation, in which VL(It=^. The other is energy dissipation 
in the inner part of the star due to differential rotation and some storage of energy in the 
convection zone. This energy could be omitted from the luminosity variation. These two 
modifications lower the luminosity variation by a factor of 4 but we are still left with a large 
value of ALijM5=3.31 L0=O.85 L. Therefore, we conclude that the Applegate mechanism 
is not sufficient to explain the observed period changes. 



8 RESULTS AND DISCUSSION 



The new 24 high-resolution echelle spectra of CF Tuc were analysed and precise spectroscopic 
orbital elements were obtained by means of two techniques; cross-correlation and spectral 
disentangling. The KOREL program was applied to four spectral orders, which contain 
about 15 lines, and then reliable radial velocities of both components of the system were 



Collier et al. 



obtain ed. I n the literatu re, there are two spectroscopic studies on CF Tuc: 
f ll98ll ) and iBalonal (119871 ). Authors of the former obtained 31 spectra of the system between 
1976 and 1978 and used absorption and Call H and K emission lines to derive radial 
velocities of both components. However, their data have a large standard error of 14 km/s, 
and t heir spe c trosco pic orbital elements were less accurate than these presented in this 



work. 



Balonal (119871 ) gave radial velocity measurements of only the hotter component and 



its orbital solution. 

Using the observed spectral lines with high precision (i.e. Na D2 line), we found the 
components of CF Tuc to be in synchronous rotation. In fact, the deriv ed rotational v elocit y 



of this secondary star is puzzling, a point which was examined by 



Coates et al 



(1200 



They noted that there is a discrepancy between measurements of the rotational velocity of 
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the se c ondary publis h ed in the hterature (IBudding and McLaughhnI (119871 ): 
(I1999I); 



Anders et aL 



DonatietaL 



(jl997l )). The rotational velocity of the primary (hotter component) 
ranges from 25 to 30 km/s, which is similar within uncertainties, while that of the secondary 
component is 52 to 70 km/s, which corres ponds to a ra t her la rge range of values in the 



secondary star radi us, from 3 Rq to 4.3 -Rq. 



Coates et al 



Donati et al. 



(120001 ) used the measurements of 



(119971 ) and emphasized that the projected rotational velocity of the secondary 
should be ~ 70 km/s from the derived absolute parameters of the system if the components 
rotated synchronously. 

Our radial velocity and the 2007 BV light curves from Innis, were simultaneously solved 
using Wilson- Devinney code. The 2007 BV light curves show large asymmetry in the two 
different maxima similar to those observed in RS CVn type eclipsing binary stars. CF Tuc 
is a bright system (V = 7.6 mag), and has been frequently observed photometrically. In the 
last three decades about 30 light curves of the system were obtained; almost all of them 
exhibit large asymmetries. The secondary component is apparently magnetically active, as 
Collier et al. (1981) observed the indicator of magnetic activity. Call H & K emission lines, 
from this component. Our Ha and Call H & K observations also show that the secondary 
component is a chromospherically (or magnetically) active star. Therefore, these light asym- 
metries we re considered as m a culati on effects and interpreted using spot models on the 



secondary. 



Budding and ZeilikI (119951 ) solved 25 light curves, mainly taken in broadband V 



and spanning a 16 year period, and modelled the spot activity of the secondary (cooler) 
component using their program (I LOT). They suggest ed that the spot luminosity has been 



decreasing over the 16 year period. 



Anders et al. 



(119991 ) solved 27 light curves of the system. 



taken between 1979 and 1996, and estimated the parameters of the spot placed on the sec- 
ondary, especially the longitudes and radius of the spots. They suggested that there was a 



strong tendency for spots to appear in a narrow range of longitudes. 



ust b efore the phase of 



primary minimum and just after secondary minimum. lAnders et al.l (11999! ) determined than 
in the 2007 BV light curves, the maculation wave begins to appear just after the phase of 
secondary minimum. Therefore, we used one dark spot on the secondary to solve these light 
curves. On the other hand, the Hipparcos light curve shows two distortion waves; the first 
one appears just after the secondary minimum like the 2007 light curves and the other is in 
the primary minima. Therefore, we introduced two dark spots on the secondary component 
to solve the Hipparcos light curve. Since the parameters obtained from the simultaneous 
solution of BV light and RV curves are more reliable, during the iterations of the Hipparcos 
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Table 8. Absolute parameters of the CAB star CF Tuc 



Parameter 


Primary 


Secondary 


M (Mo) 


l.lliO.Ol 


1.23±0.01 


R {R@) 


1.63±0.02 


3.60±0.02 


Log g (cgs) 


4.05±0.02 


3.42±0.02 


T (K) 


6100±200 


4286±219 




3.45±0.17 


3.27±0.23 


L (Lq) 


3.32±0.51 


3.91±0.84 


My 


3.50±0.17 


3.80±0.23 


Miol [system) 


2.60±0.15 


My {system) 


2.87±0.15 


d (pc) 


89±6 



light curve only spot parameters and luminosity of the primary have been adjusted. A com- 
parison of the light curve asymmetries with Ha emission modulation suggests a close spatial 
association between the photospheric spot(s) and chromospheric/coronal active region of the 
secoi idary component. Such associations have a l ready been found in several CAB system s 



[e.g. 



Frasca et al. 



fll998h : 



Catalano et al. 



Biazzo et al. 



fcood); 



Frasca et al. 



(|2008|)). 



mm -. 

The simultaneous solution of BV light and radial velocity curves allows us to calculate 
the absolute parameters of CF Tuc. The resulting parameters (with uncertainties) are given 
in TablelHl As mentioned above, following Anders et al. (1999) and Budding and McLaughlin 
(1987) we adopted an effective temperature of about 6100 K for the primary component of 
the system. Anders et al. found it to be in good agreement with the value derived from 
the colour indices of Collier (1982), and Budding and McLaughlin estimated it using the 
spectral type of the primary revised by Collier et al. (1982). If the standard error of the 
photometric observations by Collier is taken as 0.01 mag, this corresponds to an uncertainty 
in the primary component's temperature of 200 K. On the other hand, the uncertainty in 
the temperature of the secondary component, 19 K, given in Table 6, is the formal la error 
coming from the WD simultaneous solution. The corrected uncertainty could be estimated 
as 219 K based on the uncertainty of 200 K in the effective temperature of the primary. In the 
calculations, the temperature, bolometric magnitude and bolometric correction of the Sun 
were taken as 5780 K, 4.75 mag and -0.14 mag, respect ively. Bolometri c corrections for the 
components of the system were taken from the tables of IZombeckl (jl990l ). From the distance 
modulus of 2.87±0.15 mag, we derived the distance of the system to be 89 ±6 pc, under 



van Leeuwen 



the as sumption of ^4^=0. According to the new Hipparcos parallax given by 
( 120071 ) ■ the distance to CF Tuc is about 89±4 pc. This consistency between the dynamic 
and Hipparcos parallaxes shows the accuracy of the determined absolute parameters of CF 
Tuc. 
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Figure 10. Hotter and cooler components of CF Tuc together with those of the other RS CVn type echpsing binaries, in 
which the cooler component is more massive (data from Eker et al. 2008), plotted in the HR diagram (left panel). Comparison 
between evolutionary models and the physical parameters of CF Tuc in the mass-radius diagram (right panel). The filled and 
open circle symbols represent the hotter and cooler components in the binaries, respectively. The error bars of the measured 
quantities are shown by vertical and horizontal lines. Zero Age Main Sequence (ZAMS), Terminal Age Main Sequence (TAMS), 
the evolutionary tracks and isochronc were taken from Girardi et al. (2000) for the solar chemical composition. The numbers 
denote initial masses. 



The locations of the components of CF Tuc in the luminosity-effective temperature {Teff- 
L, i.e. Hertzsprung- Russell diagram) and the mass-radius (M-R) planes are shown in Fig.fTOl 
We considered only RS CVn type eclipsing binaries, in which the cooler component is a 
more massive one, to compare the CF Tuc system with other RS CVn's. According to these 
diagrams, while the cooler component has evolved behind the terminal age Main Sequence, 
the hotter one is still on the Main Sequence, approaching the TAMS. This case predicates 
that CF Tuc is similar to other RS CVn's. We compared our observationally determined 
physical parameters with those inferred from the evolutionary tracks in the HR diagram. 
The best fit appears at age of 5.012 Gyr for both components. The PARAM program at 



web page (http:/ /stev.oapd. inaf.it/~lgirardi/cgi-bin/param) based on the Bayesian method 
of da Silva et al. (2006) gives the CF Tuc age as 5.3 Gyr, which agrees well with our 
determination. 

The orbital period change of the system is interesting and difficult to solve. The O — C 
data of CF Tuc could be represented either by two abrupt period changes or by a sinusoidal 
period variation superimposed on a downward parabola. In the first case, there appears to 
be two sudden period jumps: in 1986 and 1995. The first jump shows a period increase. 
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while the second denotes a period decrease. According to the conservative mass transfer 
mode, the period increase corresponds to a mass transfer from the less massive to the more 
massive component. In the case of CF Tuc, since the less massive (primary) component is 
quite far from filling its Roche lobe, this period increase derived from the first jump can not 
be explained by the mass transfer mechanism. If we take into account the representation of 
— C sinusoidal variation superimposed on a parabola; the downward parabola indicates the 
highest rate of period decrease among the CAB systems. However, since CF Tuc is a detached 
system, we do not expect direct mass transfer between the components. Therefore, another 
possible explanation could be a mass loss by stellar wind from the subgiant component. Since 
the active, subgiant component fills ~ 89% of its Roche lobe, if a large amount of escaping 
mass by stellar wind transfers to the primary one, the remaining amount of escaping mass 
would leave the system. To check this hypothesis, the Ha and Call H & K spectra of the 
system were examined. Since the RVs of Ha and Call H & K emission features follow the 
orbital motion of the secondary (cooler) component, these emission features should come 
from chromospheric and/or coronal layers of that component. From the large width of the Ha 
emission line profile, we estimate turbulent velocities to be up to 200 kms~^. Such velocities 
well support our hypothesis that the emission features originate in the circumstellar material. 

The sinusoidal form of the orbital period variation was considered as an apparent change 
and interpreted in terms of the light-time effect due to an unseen component in the system. 
The large amplitude (~0.04 day) and small period (~18 years) of the sinusoidal form of the 
— C diagram give quite a large value (2.7 Mq) for the minimum mass of the hypothetical 
third body. The observed systemic velocity variation and Hipparcos intermediate astrometric 
data of CF Tuc partially support the hypothesis of the existence of a third body in the 
system. If such a third body were a main-sequence star, it would be a blue dwarf of a 
late-B spectral type. However, neither new high-resolution spectroscopic observations nor 
photometric analysis show any evidence to confirm the presence of such a star. Additionally, 
we could not expect such a young star as the binary system is much older. Therefore, the 
hypothetical third body must make a negligible optical contribution to the total light. With 
this mass, it could be either a massive neutron star (NS) or a black hole (BH). In the case of 
the companion b eing a compact ob j ect, th e question is, can we see it, i.e. in X-rays? CF Tuc 



was observed by 



Franciosini et al. 



(120031 ) and they detected X-ray photons with energies 
of one keV in quiescence and a few keV during the flares. This soft X-ray emission was 
identified as being released in the corona of the magnetically active component. A young NS 
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should give much higher X-ray emission , and an isolated NS would produce a termal X-ray 
emission of 40-100 keV (ITriimperl (120051 )). The progenitor of the NS must be a B-type star 
which evolves fast through the Main Sequence, within 100-500 million years, depending on 
the mass. The newborn NS would have a temperature of 10^"*^ K and could be detected 
in X-rays. However, assuming the triple system has been formed at the same time and not 
via a capture, the NS would have enough time to cool down to a much lower temperature. 
According to the NS models, after 10^ years, the NS temperature would be only 10^ K. 
The other possibility, if either a NS or a BH is the third component, could be the accretion 
luminosity when matter is being accreted onto this object. In fact, a black hole can be seen 
only through its accretion effects. Due to the large distance from the binary system to the 
companion (4.9AU), the only possibility for accretion would be through a stellar wind from 
the binary. Even though CF Tuc has the higher rate of stellar wind among CAB binaries, at 
the third body separation only a tiny fraction of mass lost by the active component could 
be accreted, as such accretion is much less efficient than tha t through the Lagr angian point. 
We shall estimate the rate according to the formula given by lFrank et al.l (120021 ). Taking 6M 
= 3.38 X 10"^ Mq/jt, separation of 4.93 AU, mass of the accreting object to be 3 Mq and 
derived radius and mass of the secondary star, we calculated the accretion luminosity to be 
about Lace = 3.4 X 10^4 erg/s. Suc h a va lue is comparable to the quiescent luminosities of 
NSXN and BH XN (INarayan et al.l (120021 )) and quiescent low-mass X-ray binary transients 
( Lasota ( 200ol )). ^f the third companion is a NS with such a luminosity it should be observed 
in X-rays, unless the efficiency of the mass transfer to the primary is more efficient than we 
have assumed. If the companion is a BH, which has no surface, a significant fraction of mass 
could be lost below the event horizon, instead of being converted into hard photons. Also, 
if the accretion onto the BH is spherical and not through a bow shock, as considered above, 
it could be less efficient as well. How ever, the fac t ors de s cribing this effi ci ency are withi n 
a huge range betw een 10"^ to 10"^ ( Frank et al. (2002), ShapirJ (1974). Shapiro ( 1973 ). 



Petrich at al. 



(119891 )) and strongly model dependent (i.e. rotation of the BH, magnetic field, 
speed of the BH or the accreted material). Concerning the efficiency of the accreting matter 
being converted into photons, this factor could be in the same range. If the factors are in the 
lower end, then the existence of a BH as the third companion in CF Tuc in not impossible, 
though unlikely, as the possibility of its observation may fall below the detection limit and 
the light time effect would be the only evidence of its existence. 

Another explanation of the orbital period sinusoidal variation could be the Applegate 
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mechanism. However, the Applegate model with Ms=0.1M, fi^^ = and ALjims = O.IL 
cannot explain the rate of orbital period change observed in CF Tuc. 

The last possibility, we can consider, would be a spurious sine term which may not repeat 
in the future. The observed light curves have strong asymmetries which can affect accuracy 
of the minima times determination. In some cases, these asymmetries are so large that sec- 
ondary minima can not even be easily recognized. Therefore, so me spurious O — C re s iduals 



Budding and ZeilikI fll995h 



may come from the asymmetries of the light curve. For instance, 
found an approximate 6 year magnetic cycle, considering that these light curve asymmetries 
were caused by maculation effects. Therefore, we conclude that only future observations can 
reveal the true nature of the observed period changes. If they are caused by a companion 
BH, this will be the closest black hole to the Earth. 
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